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Abstract 

We show that the one-dimensional Oslo rice-pile model is a special 
case of the abelian distributed processors model. The exact steady 
state of the model is determined. We show that the time evolution 
operator W for the system satisfies the equation W n+1 = W n where 
n = L(L + l)/2 for a pile with L sites. This is used to prove that W 
has only one eigenvalue 1 corresponding to the steady state, and all 
other eigenvalues are exactly zero. Also, all connected time-dependent 
correlation functions in the steady state of the pile are exactly zero for 
time difference greater that n. Generalization to other abelian critical 
height models where the critical thresholds are randomly reset after 
each toppling is briefly discussed. 



1 Introduction 

In their pioneering work in 1987, Bak et al proposed the well-known sandpile 
model as a paradigm for self-organized criticality pQ. Since then, many dif- 
ferent variants of sandpiles have been studied. These include models where 
the height of the pile is a real variable [2], different toppling rules [3], with 
preferred direction stochastic topplings 0, models with fixed energy [H], 
etc.. While there is a fair amount of understanding by now about many of 
these models, the understanding of different universality classes of critical 
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behavior possible, and of relationship to the critical steady states seen in 
other driven systems remains unsatisfactory [7j. 

The static properties of the BTW sandpile model can be related to the 
equilibrium properties of q — > limit of the g-state Potts model and the 
statistics of spanning trees jH] . The directed model [I] can be related to the 
well-known voter and the Takayasu aggregation models However, the 
exponents characterizing the avalanche distributions in sandpile models are 
not determined completely by the critical exponents of the q = Potts model. 
We have argued that the generic behavior of sandpile like models is in the 
universality class of directed percolation JU]. However, the critical behavior 
of special models, like the BTW or the Manna models, is certainly not in this 
class. Paczuski and Boettcher have given arguments relating the exponents 
of sandpile model to that of growing interfaces ^1] . Generically, the growing 
surfaces are expected to be in the universality class of Kardar-Parisi-Zhang 
type models, not of directed percolation. It seems desirable to have a closer 
look at the relation between sandpile models and growing surfaces. This 
relationship is most clearly seen ^2] i n the simple Oslo rice-pile model 

The Oslo rice-pile is one of the simplest of models of self-organized crit- 
icality and seems to be able to describe quite well the real experiments on 
rice-piles |14j . It shows non-trivial avalanche exponents even in one dimen- 
sion. There are some exact results known about this model ^3], and the 
values of avalanche exponents are known quite accurately jTHj. The damage 
propagation in the model has also been studied |T7j. It seems worthwhile to 
try to see if the model can be solved exactly. While attempts in this direc- 
tion have not been successful so far, it does seem to have a rather unusual 
mathematical structure. In this paper, I will try to outline some of these 
properties, and hope that this will encourage further work. 

2 Definition of the Model 

The Oslo rice-pile model is a one-dimensional cellular automaton model with 
stochastic toppling rules. It is defined as follows ^3]: We consider a line of 
L sites, labelled by integers 1 to L. At each site i, there is a non-negative 
integer height variable hi, called the height of the pile at that site. We start 
with a configuration in which all heights are zero. At each site i, there is 
also a variable af, called the critical threshold, which takes values 2 or 1 with 
probabilities p and q = (1 — p) respectively, independently at each site. 
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We shall call hi — hi + \ as the local slope of the pile at site i, and denote it 
by <7j. We define h^+\ = 0, so that the local slope at the site i = L is equal 
to hi- We call a site % stable if <7j < of. A configuration is said to be stable, 
if all sites in the configuration are stable. 

If a site i is unstable, it relaxes by sending one grain to site on its right. 
In this process, hi will decrease by 1, and h i+1 will increase by 1. If a 
configuration has more than one unstable sites, then they are updated in 
parallel. ( We will show later that actually the order of topplings does not 
matter). After each toppling at a site, the critical threshold at that site 
is reset, and takes a new random value 1 or 2 with probabilities q and p 
respectively. 

The system is driven by adding a grain at the leftmost site % — 1. This 
increases the local height by 1. If this leads to an unstable configuration, it 
is relaxed by topplings till a stable configuration is reached. And then we 
again add a grain at i — 1. 

If we evolve the system like this, it reaches a self-organized critical state 
at large times, where the distribution of avalanches shows a power law tail 
with an L-dependent cutoff. 

3 The Abelian Property 

It is convenient to specify the configurations in terms of local slopes Oj. A 
particular configuration {hi} corresponds to a unique set of slope variables 
{<7;}, and vice versa. In terms of the slope variables, the evolution rules of 
the rice-pile become the following: 

i) If a particle is added from the left, U\ increases by 1. 

ii) a. If Oi > of with 2 < i < L — 1, then Oj decreases by 2, and Oj_i and a i+1 
increase by 1. 

b. If o\ > of, then, o\ decreases by 2, and o"2 increases by 1. 

c. If ctl > a c Ll then, o~l decreases by 1, <tl-i increases by 1. 

iii) After a toppling at site i, of is reset to a new randomly chosen value. 

It is interesting that these rules are same as for a one-dimensional critical 
height model, where the critical threshold is a random variable, and its value 
at a site is reset randomly after each toppling at that site. This correspon- 
dence to an abelian critical height model is not possible in higher dimensions, 
or if grains in the Oslo model are added at all sites, and not only at the left 
end. 



3 



It is convenient to work with this one-dimensional critical height model, 
in which the "particles" are the "slopes variables" of the original model. To 
avoid confusion, in this paper we shall call the rice-particles of the original 
model "grains" , and use the word "particles" to refer to the slope variables. 
Thus, in a single toppling event at site i, (2 < i < L— 1), a grain is transferred 
from i to (i + 1), but two particles are moved from i, one each to the sites 
(i — 1) and (i + 1). The topplings at the boundary sites can be similarly 
described. 

As another simplification, we need not keep track of the critical thresh- 
olds, if we adopt the following rule for relaxation whenever a new particle is 
added to a site : If, after addition, the new height is > 2, the site topples; if 
the height is 2, it topples only with probability (1 — p); and if the new height 
is 1, it does not topple. Clearly, the actual evolution is same as under the 
original definition. 

This model, then is a special case of the general Abelian distributed 
processors model defined earlier 0. We consider each site as a finite-state 
automaton, with a local pseudo random number generator (PRNG), and a 
single integer giving the number of particles at the site. If an added particle 
makes the number of particles greater than the current threshold, a fixed 
number of particles, as specified by the matrix A, are transferred to other 
sites. A new random number is drawn from the PRNG, and used to choose 
a new threshold. It is easily shown that in any configuration with more 
than one unstable sites, the final state of this system does not depend on 
the order in which these sites are relaxed, for any particular sequence of 
thresholds generated by the PRNG. This establishes the abelian property for 
the rice-pile model. 

4 Operator Algebra 

Let us denote by V a linear vector space spanned by its basis vectors \C), 
where \C) are the stable configurations of the rice-pile. We shall characterize 
the configurations C by the slope variables There are 3 L such configu- 

rations for a chain of length L. Note that the values of the critical slopes of 
are not specified in C. 

The state of the rice-pile at time t is characterized by a probability vector 
\Prob(t)) , which is an element of V. If the probability that the rice-pile occurs 
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in configuration C at time t is given by p(C, t), we write 

\Prob(t)) = Y,p(C,t)\C). (1) 

c 

We now define linear operators ( i — 1 to L ) acting on V, by their 
actions on the basis vectors \C). Starting from the stable configuration C, 
we increase the slope variable cr, by 1, and allow the resulting configuration 
to relax if necessary. If the resulting configuration is C with probability 
Wi(C'\C), we define 

a t \C) = ^(C'\C)\C). (2) 

a 

We then have 

[cii,aj] = 0, for all (3) 

The time-evolution of the system is Markovian. Let \Prob(t)) be the proba- 
bility vector of the state of the system after t paticles have been added, and 
the system is allowed to relax. We then have the master equation 

\Prob(t + 1)) = W \Prob(t)), (4) 

which defines the time evolution operator W. For the rice-pile model, with 
particles added only at the left end, clearly, we have W = a\. 

Let us now give a more explicit representation of the operators {d}. If 
Oj acts on a configuration where 0; L = 0, we get a state with <jj increased to 
value 1, and all other cr's are unchanged. 

o i |...,<T i = 0,...) = | . . . ,(Ti = 1, . . .). (5) 

If a* acts on a configuration with — 1, then with probability p, it increases 
its value to 2. With probability (1 — p), it causes a toppling there, which 
would change <j{ to 0, and add a particle at sites (i ± 1). 

Oi|...,(7i = 1,...) = p| . . . ,0i = 2, . . .) + qa i - 1 a i+1 \ . . . , a* = 0, . . .). (6) 

Acting on a configuration with Oi = 2, d will always cause a toppling, and 
we have 

Oj| ■ ■ ■ ,<7j = 2, . . .) = dj-iOj+il . . . ,<7j = 1, . . .). (7) 

These equations fully define the action of the operators aj. They also 
hold for the boundary sites i — 1, and i = L, if we assume the conventions 

a = 1, and a L+i = a L . (8) 
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Applying these rules repeatedly, we can determine the effect of any of the 
the operators {a^} on any stable configuration. For example, consider the 
case L = 2. Let us determine the result of a\ acting on |2, 2). We get, 

ai|2,2) = a 2 |l,2) = aia 2 |l,l) = a x { 2) + q ai \l, 1)). (9) 

These two terms can be evaluated further as 

ai|l,2)=p|2,2) + ga 2 |0,2), (10) 

with 

a 2 |0,2) = aiO2|0,l) = a 2 |l, 1) = p\l, 2) + gai|l, 1), (11) 

and 

ai|l, 1) = p|2, 1) + ga 2 |0, 1) = p\2, 1) + pq\0, 2) + g 2 |l, 1). (12) 
Putting all these together, we get 

ai |2, 2) = p 2 |2, 2) + (1 - p 2 ) [p|l,2> + gp|2, 1) + q 2 p\0, 2) + g 3 |l, 1)] (13) 

We also note that these operators satisfy the equations 

af = aj_iaj<2j +1 , for all i. (14) 

This is equivalent to the observation that on adding three particles at site i, 
a toppling there must occur at site i, so that it is same as adding a particle 
each at the sites i — 1 , i and i + 1 . 

Given any product of the form 02^1070202 . . ., we first use the commuta- 
tivity property to bring it to the form a^a^ 2 . . ., where ni,ri2, . . . are non- 
negative integers. Then, we use the reduction rules Eq. (JUjl to express it 
as a sum of lower powers such that n, < 2 for all i. Consider evaluating a\ 
for r = 1, 2, . . .. As the number of possible answers is finite, we must have a 
minimum value of r, such that a\ is equal to a lower power of a±. Normally, 
one would guess that r is of order exp(L), the number of different terms al- 
lowed. Interestingly, it turns out to be much smaller, and we now show that 
r = L(L + 1) + 1 for all L. 

Using Eq.((T2J), it is straightforward to show that 

L(L+1) 2 2 2 2 /-I r\ 

a 1 — a 1 a 2 a 3 . . .a L . (15) 

Multiplying both sides of this equation by ai, and again using the reduction 
rules Eq.()14p. we get 

L(L+1)+1 2 2 2 2 

% =a x a 2 a z ...a L . (16) 
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As ai = W, we get 

yyi(L+l) + l = y\/L(L+l) ^ 

This equation holds as an operator equation, over the 3 L dimensional 
vector space V, the space of all stable configurations. This has the remarkable 
consequence that eigenvalues of W are either or 1. W is a stochastic 
matrix. As there is a finite probability that the maximal state |2, 2, . . . , 2) 
can be reached from any stable configuration, all recurrent configurations are 
reachable from any other, and there is a unique steady state for this Markov 
process. This implies that there is only one eigenvector of W with eigenvalue 
1. Let us call it \^ s t) Thus all the other 3 L — 1 eigenvalues of W are 0. 

From Eq.(|17|). it follows that for any vector \<p) in V, W L ^ L+1 ^\(f)) is pro- 
portional to \^st)- Also, as W preserves the probability sum, for any initial 
basis vector \C), we get 

A simpler way to determine \^ s t) is to use the equation 

|* st ) =^2,2, 2... 2). (19) 

This equation is interesting, as it says that if we take the initial config- 
uration as the one in which each Oi has the highest allowed value, and just 
add one more grain, the stochastic evolution would result in different final 
stable configurations with probabilities exactly equal to their values in the 
steady state. To prove this result, we need only note that for any stable 
configuration, 

a i \...,a i = 2,...)=a 3 i \...,a i = 0,...). (20) 
Then we can write 

Oi|2, 2, 2 ... 2) = al\0, 2, 2 . . .) = aia 2 \0, 2, 2 . . .) = a^o, 0, 2 . . .) (21) 
= ala 2 a 3 \0, 0, 2 . . .) = a\a 2 a\\Q, 0, . . .) = etc., (22) 

finally giving 

ai|2,2,2...) = alala 2 3 ...a 2 L \0,0,0...) = \^ st ). (23) 

In fact, it is possible to get a result stronger than Eq. (fTH|) . We note that 
most of the stable configurations are transient, and do not occur in the steady 
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state. The different recurrent configurations of the Oslo rice-pile have been 
characterized ^H]- For any recurrent configuration C with m grains, arguing 
as above, we can show that 

a? L+1) ~ m \C) = a\a\a\ . . . a||0, 0, . . . 0) (24) 

Since the lowest allowed value of m in recurrent configurations is L(L + l)/2, 
corresponding to the configuration |1, 1, 1, . . .), we get that for any recurrent 
configuration 71, 

4 (L+1)/2 |^> = I**) (25) 

Let X(t) be some scalar observable whose value depends on the stable 
configuration of the pile after t grains have been added. Then any particular 
evolution of the rice-pile will generate a stochastic time-series {X(j)},j = 
1,2,.. .. For example, X(t) may be the total height of the pile h\ at time t, 
or the number of grains in the pile. Clearly, this time-series has nontrivial 
correlations. For example, hi(t+l) < hi(t) + l. For any two such observables 
X(t) and Y(t), we define the connected time-dependent correlation function 
Cxy(t) in the steady state as 

C XY (r) = (X(t)Y(t + r)> - (X(t))(Y(t + r)> (26) 

However, from Eq.()25|>. it follows that the conditional expectation value of 
Y{t + r), given that the observable X has a particular value X(t), must 
be equal to its unconditional expectation value, whenever r > L(L + l)/2. 
Thus, Cxy(j) is exactly equal to zero for r > L(L + l)/2, for all observables 
X and Y. 

For t < L(L + l)/2, the correlation function Cxy(t) is not always zero. 
This correlation function can be nontrivial, even though W has only one 
nonzero eigenvalue, because W is non-hermit ian. This is easily checked for 
the simple cases L = 2, 3. 



5 Generalizations 

It is straightforward to generalize the discussion given above to more general 
abelian sandpile model (ASM) Consider a stochastic sandpile model 

defined on a set on iV vertices labelled by integers 1 to N. At each site 
is a height variable Zi, which takes non-negative integer values. There is 
an iV x Af integer matrix A, which specifies how particles are transferred 
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under topplings: If Zj is greater than or equal to a site-dependent threshold 
value Zj C , there is a toppling at site i, the heights at all sites j are updated 
according to the rule 

Zj -> Zj - Aij. (27) 

The matrix A is assumed to satisfy the good-behavior conditions, same as for 
ASM [10]. After each toppling at a site i, the critical threshold at that site is 
randomly reset to a new value Zi C = r, with r lying in a finite range of values 
z™ n and from a known probability distribution Probi(zi C = r). This 
distribution functions can be different at different sites. Note that while the 
threshold is randomly reset, we are assuming that the matrix A specifying 
how particles are transferred under toppling does not change. Without any 
loss of generality, we may assume that z™ m = A^. This corresponds to 
setting the minimum allowed value of Zi to be zero for all sites i. 

It is easy to see that this model is still abelian, and that just like the 
ASM with deterministic rules, there are forbidden sub configurations in this 
model. In fact the forbidden subconfigurations are exactly the same as in 
the deterministic ASM specified by the toppling martix A [TSJ. The number 
of allowed configurations in the stochastic model is larger, as the maximum 
allowed value of Z{ is larger. We define a new matrix A such that 

&ij = Aij, for i jLj- (28) 

and 

hi = <T • (29) 
Then, it is easy to see that any allowed configuration for a deterministic ASM 
with toppling matrix A is also a recurrent configuration of the stochastic 
model, and vice versa. The number of recurrent configurations is given by 
det(A). However, all configurations are not equally likely to occur in the 
steady state. 

For the 1-dimensional model studied in this paper, the only non-zero 
entries of A are A^j = 3, for i < L, and A^ ^ = 2, and A iA ±i = —1. The 
determinant is easy to calculate, recovering the results of [T5] . 

Again, for this stochastic model, we define operators {a{\ corresponding 
to adding a particle at site i, and relaxing. Then these operators commute 
with each other. In addition, they satisfy the equations 

., max ( ^max ^.min \ * 

a>i = a i ' YY a j ■ (3°) 
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This equation is a generalization of the Eq. (jl4j) in the previous section. 
As before, we find that the eigenvalues of the operators {cij} are either all 
zero, or same as that for the deterministic ASM with toppling matrix A. 
The steady state \^f s t) of the stochastic model is given by 

/ -max ~min\ 

l^> = n«r' c ~ lc W, (si) 

i 

where |$ s t) is the steady state of the deterministic ASM with toppling matrix 
A. To prove this statement, we note that if for the deterministic model 
a>i\C) = \C), then for the stochastic model we must have 

( max y min \ I -max „min\ 

aiU a i \C) = U a i (32) 

i i 

as each toppling which is required to relax from C to C in the deterministic 
case, will also be allowed in the stochastic case with the larger threshold 
values, so long as there are enough extra untoppled particles present [TS]. 
Then, if |<J> st ) is an eigenvector of all aj with eigenvalue 1 in the deterministic 
model, \^st) would be so for the stochastic one. 

It may be noted that this operator algebra does not depend on the proba- 
bility distribution of the random thresholds, and depends only on the values 
of z?]? x and z™ n , and the toppling matrix A. All nontrivial eigenvalues of 
{a{\ are solutions of algebraic equations (|30|) . and are the same as for the 
corresponding ASM with toppling matrix A. 

The Oslo rice-pile model is special in that in this case, the deterministic 
ASM with same matrix A has only one recurrent configuration, and the 
= |1, 1, 1, . . .). This explains why W has only one non-trivial eigenvalue. 

The distribution of avalanche sizes, and the weights of different configura- 
tions in the steady state certainly do depend on the probability distribution 
of thresholds. These are rather difficult to calculate explicitly even for the 
original one-dimensional model, except for very small values of L. Exact 
calculation of these for general L seems more difficult. This seems to be an 
interesting open problem. 
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